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Abstract: The Hanbury-Brown Twiss correlation function for two identical parti- 
cles is studied for systems with cylindrical symmetry. Its shape for small values of 
the relative momentum is derived in a model independent way. In addition to the 
usual quadratic "side", "out" and "longitudinal" terms in the exponent of the cor- 
relator, a previously neglected "out-longitudinal" cross term is found and discussed. 
The model-independent expressions for the size parameters of the HBT correlation 
function are interpreted as lengths of homogeneity of the source, in distinction to 
its purely geometrical size. They are evaluated analytically and numerically for two 
specific thermal models featuring collective transverse and longitudinal flow. The an- 
alytic expressions derived allow one to establish qualitatively important connections 
between the space-time features of the source and the shape of the correlation func- 
tion. New ways of parametrizing the correlation function and a new approach to the 
measurement of the duration of the emission process are suggested. 
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1 Introduction 



It is widely accepted that if the nuclear matter created in ultra-relativistic heavy-ion 
collisions attains a high enough energy density, it will undergo a phase transition 
into a quark-gluon plasma. For this reason, it is of great interest to determine the 
energy densities actually attained in these collisions. The total interaction energy of 
a given reaction can be directly measured by particle calorimeters and spectrometers. 
Although there is no analogous direct measurement for the size of the reaction region, 
Hanbury-Brown Twiss (HBT) interferometry [|TJ provides an indirect measurement for 
both the spatial and temporal extent of the reaction region in terms of the correlations 
between produced particles. 

Consequently, the greatest challenge for theorists studying HBT interferometry 
today is to determine exactly what information the reported experimental correlation 
radii are telling us about the source. Obviously, the most powerful statements to 
this effect are those which can be made in a model-independent fashion. Although 
the individual reactions measured experimentally may not be completely cylindri- 
cally symmetric, it is safe to assume that a large ensemble of similar reactions will 
produce cylindrically symmetric data. For this reason, we have generalized the work 
of f|| by using the covariant Wigner function formulation [2-6] of HBT interferome- 
try to derive cylindrically symmetric, but otherwise model independent expressions 
for the correlation radii, both using standard cartesian momentum differences and 
boost-invariant rapidity differences. Two important model-independent statements 
can then be made. First, cylindrical symmetry in no way precludes the existence of 
an "out-longitudinal" cross term in the correlation function [/J, and in fact in general 
such a term would be expected to appear. Second, the correlation radii do not nec- 
essarily measure the geometrical size of the reaction region, but rather the lengths of 
homogeneity of the source as seen by a particle emitted with the average momentum 
of the studied pair || . 

To see how these effects manifest themselves in a concrete (though still qualitative) 
way, we apply our model-independent formalism to two specific thermal models, both 
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of which feature a constant freezeout temperature. The first model is a generalization 
of ||, featuring nonrelativistic hydrodynamic flow which, however, can be different 
in the longitudinal and the transversal directions. Since this model is completely 
gaussian, it is easy to verify explicitly that the spatial lengths of homogeneity depend 
not only on the geometrical size of the reaction region, but also on the spatial gradients 
of the hydrodynamic flow. Similarly, the cross term just measures the temporal length 
of homogeneity, which in this nonrelativistic case is simply the duration of particle 
emission. 

The second model that we consider is a variation of [|l(J, featuring a Bjorken 
scaling longitudinal flow and a nonrelativistic transverse flow. Although this model 
is not completely gaussian, analytic results derived from a modified saddle point 
approximation are able to reproduce numerically generated results to within 20-30% 
for pions and much better for kaons. The analytic results provide valuable qualitative 
insights into the generic influence of various physically relevant parameters of the 
source distribution on the shape of the correlation function. We show that this model 
features a large cross term whose effects can clearly be seen in a two-dimensional 
plot of the "out-longitudinal" correlation function. In addition, we show that the 
theoretical interpretation of the correlation radii simplifies immensely when rapidity 
differences rather than longitudinal momentum differences are used to parametrize 
the correlation functions. In light of these results, we make explicit suggestions of 
useful new ways in which experimentalists can organize their measured correlation 
data. 

2 Model Independent Correlation Radii 

The HBT correlation function for two identical on-shell particles is given by 11 



mass m and 3-momentum p, P2 is the corresponding invariant 2-particle distribu- 





where P\{p) = E p (dN/d 3 p) is the invariant 1-particle distribution for a particle with 
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tion function, and N {N 2 ) is the average number of particles (squared) produced 
in a reaction. By quite general arguments it can be shown that in the plane wave 
approximation for chaotic sources [2-6] 



Jd 4 xS[x i(pi +p 2 )}e iq - x 2 

Pl{Pl)Pl{P2) 

where the + (— ) sign is for bosons (fermions), q = pi — P2 is the 4-momentum 
difference of the two particles, and p\ = Ei are the on-shell energies. Furthermore, 
the emission function S(x,p) is a scalar function of the 4- vectors x and p which obeys 

Jd A xS(x,p l ) = P 1 (p t ) . (3) 

As an example, in the local hydrodynamic formulation involving a sharp 3-dimen- 
sional freeze-out hypersurface one has jl2) 



„, x _ 1 p-n(x) 

[X,P) ~ (2tt) 3 exp[/3(x) (p-u{x) - fi{x))] T 1 ' 



where /3(x), fi(x) and 

n M (x) = Jj 3 a,(x')5^(x-x') (5) 

denote the local hydrodynamic flow velocity, inverse temperature, chemical potential, 
and normal-pointing freeze-out hypersurface element, respectively. 

2.1 Cartesian Momentum Coordinates 

In order to simplify computation, the correlation function is often approximated by 



using on-shell momenta in the emission function [fj, [|, [12], |H| . For example, one can 
define § 



where = +P2) and Kq = Ek = ym 2 + \K\ 2 . Neither the present definition of 
K nor the different definition we will use in the next subsection should be confused 
with the usual off-shell definition of K = \{Ei + E 2 ) which is suggested by eq. (0). 



We begin by using the conventional HBT cartesian coordinate system which is 
defined as follows: The "longitudinal" or z (subscript L) direction is defined to be 
parallel to the beam; the "out" or x (subscript _L) direction is parallel to the com- 
ponent of K which is perpendicular to the beam; and the "side" or y (subscript s) 
direction is the remaining transverse direction. For \(j\/Ek < 1, we then have 

q ■ x ~ (3-qt — q±p cos — q s p sin — q L z , (7) 

where p = y/x 2 + y 2 , <fi = tan~ 1 (y/x) and (3 = K/Ek is the velocity of a particle with 
momentum K. 

To present their data, experimentalists use these coordinates in one of two different 
reference frames, both of which can be obtained by a longitudinal boost from the lab 
frame: The fixed observer frame is usually taken as the rest frame of the participant 



center of mass and is the same for all particle pairs |L4|, [15|, |16[. The "LCMS" 
(longitudinally co-moving system) frame, on the other hand, is defined as the frame 
in which Kl = and thus varies for pairs with different longitudinal momentum 



in the fixed observer frame |T7| |. Consequently, as pointed out in ]TS 



TUfl , a qL- correlation function should then only be measured at a given value of K^, 
and an averaging over Kl should be avoided. However, since different values of the 
longitudinal component of the mean momentum lead to different reference frames, 
the interpretation of a possible -ft^-dependence of the correlation radii turns out to 
be conceptually nontrivial in the LCMS. Later, however, we will show that for the 
special case of a system which is undergoing Bjorken longitudinal expansion, the 
LCMS radii are nothing more than approximations of fixed frame radii which are 
evaluated in rapidity coordinates (see next subsection). To avoid the complication of 
shifting reference frames, we perform all of our calculations in a fixed frame, though 
we do point out how to find the LCMS results. 

Due to the symmetry C(pi,p2) = C(p2,pi) and the fact that when q — > the 
correlation function C(pi,p2) — *■ 1 ± 1 (as can be seen from eq. (Q)), it is reasonable 
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to assume that for sufficiently small momentum differences q, C takes the form 



c(pi, p 2 ) = f ± ex P - J2 £R\ - 2 E m^l 



(8) 



i i^j 



where the coefficients R\ and R\j depend on the average pair momentum if = + 



f>2). Note that the i?f are always positive, but the i?? can be either positive or 
negative; we simply use the it!? notation to denote the fact that they are coefficients of 
terms which are quadratic in q^. Furthermore, in order for the peak of the correlation 
function to be located at q — 0, it must be true that for all i and j 



Below we obtain model independent expressions for the "radii" R\ and i??- by effec- 
tively taking second derivatives of the correlation function with respect to and qj 
around q — 0. 

Before proceeding, we would like to point out that one must take care when 
comparing the above radii to experimentally measured correlation radii since the 
former measure second derivatives of the correlation function around q — 0, while 
the latter are parameters of a gaussian fit to the whole correlation function [14-17] 
and are essentially determined by its width. Nevertheless, there are many interesting 
"gaussian" models for which the two different ways of defining the radii give roughly 
the same results. To the extent that the part of the correlation function measured 
by experimentalists is roughly gaussian, certain of these "gaussian" models should be 
able to provide good descriptions of the data. In this work we are therefore restricting 
the application of our model independent results to "gaussian" models for which the 
simple expressions that we generate below provide valuable insights as to how various 
parameters of a given source distribution will qualitatively affect measurable features 
of the correlation function. 

Since S(x,p) transforms as a scalar under Lorentz transformations, it can be taken 



2\R ij \ 2 < R^ + R] . 



(9) 
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to have the following functional form 

S(x,p) = S(x, p-p, p-u(x), p-v(x), p-w(x), ...) , (10) 

where u, v, w, etc. are space-time dependent local 4-vectors. Cylindrical symmetry 
can be enforced by demanding that S has no explicit dependence and that all of 
the relevant local 4-vectors be cylindrically symmetric. For example, 

u{x) = (uq, u p cos0, u p sin 0, u z ) (11) 

where Uq, u p and u z are all independent of 0. Given these definitions, 

K-u = EkUq — K±u p cos — Klu z , (12) 

so for cylindrically symmetric systems S(x, K) is even in <p. 

Using (0), we can now expand the factor exp(iq-x) in eqn. (EH for small q, keeping 
only the terms even in 0, because the odd terms vanish upon integration. We find 

C(q,K) = l±{l - q 2 s {y*)-([q ± ( x -p ± t)+q L (z-l3 L t)} 2 ) 

+ (q±(x-(3 ± t)+q L (z-(3 L t)) 2 + 0[((q-x) A )}}, (13) 

where x = p cos 0, y = p sin 0, and we have introduced the notation 

(0 = ^-^Jd 4 xZS(x,K). (14) 

Eq. QT^D generalizes similar results obtained in for a 1-dimensional situation. 
Exponentiating (|13D , we can see that for any cylindrically symmetric system the 
correlation function for small momentum differences will take the form 

C(q, K) ~ 1 ± exp [-q 2 s R 2 s - q 2 ± R 2 ± - q 2 L R 2 L - 2q L q L R 2 ±L \ . (15) 

The R 2 which correspond to the approximation (||) can simply be read off as the 



coefficients of the corresponding q^qj terms in eqn. (13 ) 

?2 I „,2 



K = {v 1 

R\ = ((x - (3±t) 2 ) -(x- (3 ± t) 2 

Rl = ({z - M 2 ) - (z - M 2 

R\ L = ({x-/3 ± t){z-/3 L t))-(x-P ± t)(z-/3 L t) . (16) 



They are functions of K due to the ^-dependence of S(x, K) in definition (fH|) of the 



expectation value (...). One of the most interesting features of (|T|) is, as pointed out 
in 0], the occurrence of a q±qL cross term which has never before been discussed in 
the literature. 

Before exploring the implications of this term, we would like to give an intuitive 
interpretation of the model-independent expressions (ITW). To this end we follow the 
work of H and introduce the concept of a length of homogeneity. We begin by 
defining the spacetime saddle point x of the emission function S(x, K) through the 
four equations 



JL 

dx. 



In S{x, K) 







(17) 



where fi = {0, 1, 2, 3}. Essentially the saddle point is that point in space-time which 
has the maximum probability of emitting a particle with momentum K. A saddle 
point approximation for S(x, K) can then be made in the following way 



S(x, K) ~ S(x, K) exp 



2X l 



il>V 



where we define the length of homogeneity of the source in the /zth direction by 



and 



d 2 



dx 1 ., 



\nS(x,K) 



-1/2 



d d 
^ doc y 



\nS(x,K) 



(19) 



(20) 



From flT|), it can be seen that the length of homogeneity provides a measure of 
the region over which the source is relatively constant as seen by a particle with 
momentum K. Obviously, if a source has large temperature or flow gradients, the 
length of homogeneity may be determined by these more than by geometrical density 
gradients. 
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Notice that if B^ v C (which can always be arranged by making the right 

choice of variables), then 

K) - (x,) 2 * XI (21) 

where we do not use the summation convention. Since all of the radii of eq. (p~6|) 
contain terms of the above form, these radii are evidently measuring lengths of homo- 
geneity rather than strictly geometrical sizes. For example, the "side" radius measures 
A2 (in the y direction), which by cylindrical symmetry must be equal to the length 
of homogeneity in the transverse or radial direction. Later, we will see how these 
lengths manifest themselves in definite models. 

Corrections to the radii of (O) can be calculated by considering the exact corre- 
lation function (^) rather than the approximation (^). (Within their one-dimensional 
model these corrections were also found in 0.) The corrections to the denominator 
can be found by noticing that due to cylindrical symmetry, P\(p) is really only a 
function of the longitudinal and radial components of p. Hence, 

P\{Pi) = h\ViL, \JprPi 

= Pi (k l + \q L , K ± + \q± + J- + OQq \*/E 2 K )j . (22) 
Keeping only up to quadratic corrections in q, 



PiiPi 



d i q L d , I ( q± d ^q L d 



If q 2 s 

+ 2 { q± + AK ± I dK ± ' 2 dK L ' 2 V 2 dK ± ' 2 dK L 



Pi{K) 



(23) 

Pi{f>2) can be found simply by letting q — > — q in the above expression. Combining 
these, re-exponentiating, and again keeping only terms up to second order in we 
get 



A (Pi) Pi{n) * [Pi{K)Y exp -qj 6Rj - q 2 ± 5R 2 ± - q\ 5R 2 L - 2q ± q L 6R 2 ±L , (24) 
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where 



8Rl 



1 d 



lnPi(if) 



k l 



5ft 



1 d 2 
AdlCi 



In Pi (if) 



K L 



5ft 2 



1 d 2 
AdRj 



In Pi (if) 



5ft 2 



1 d 



4dif L Idifj 



In Pi (if) 



(25) 



Notice that all of these corrections are direct experimental observables. For example, 
5R\ is the curvature of a plot of lnPi(p) as a function of p± for fixed pi- 

Finally, we turn to the corrections induced by using the correct off-shell energy 



|(Pi +P2) i n the emission function of eqn. (0) rather than the approximate on-shell 
value Ek of eqn. @. Again making a Taylor expansion for \q | <C Ek, 



E 



K 



77T 



|if| 2 + |(|g| 2 



we can see that 



m 2 + \(\q\ 2 -(f3-q) 2 



(26) 



(27) 



Therefore, we can expand around the on-shell momentum if in the following way: 

d 



S (x, |(pi + p 2 )) ~ |l + i (|g | 2 - 0-q ) 2 ) 
To quadratic order in q, then 

' d i xS(x,\(p l +p 2 )) e iq - x ~ 



dm 2 



S(x,K) 



28 



exp 



id 



dm- 



In Pi (if) 



d 4 x S(x, if) e 



iq-x 



(29) 
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Putting everything together, we find the following corrected model-independent 
expressions for the correlation radii of eqn. (|T5|): 



+ (iA-5 (1 -^) b, ^ ; 



«i t ~ ((x-P ± t)(z-P L t))-(x-/S ± t)(z-P L t) 

+ (\d£m +1 2 p ^)^ (g) - (30) 

Note that LCMS radii can be found from the above expressions simply by setting 

0L = O. 

The first thing to observe about the above radii is that cylindrical symmetry alone 
does not cause R\ L to vanish, so a q^qi cross term (as in (|T5j)) should be included in 
any experimental fit to the data. However, it is interesting to note that for the case 
K± = (/3± = 0), S(x,K) is independent of 0, so R\ L does vanish (see appendix). 
Furthermore for this case R 2 S = R\ as it must, since if K± = it is impossible to 
define a difference between the "out" and "side" directions. This means that the q±qL 
cross term (as well as the difference between R\ and R 2 S ) will be most noticeable for 
pairs with large K±. We would also like to point out that the cross term vanishes for 
spherically symmetric systems if one redefines the z direction in the direction of K 
0, since in this case K± = by definition. For any collision experiment, however, it 
is best not to make this redefinition, since only cylindrical symmetry about the beam 
can be assumed. It should also be noted that if future heavy ion experiments are able 
to generate HBT correlation functions from a single event, then cross terms involving 
q s q± and q s qi should be included in any fits as tests of the cylindrical symmetry of 
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the individual reaction under consideration. 

Before going on, we would like to say a few words about the validity of the approx- 
imation of eqn. (||) and the size of the correction terms. Since the ffif of fl25|) can be 
measured from single particle distributions, a model-independent experimental esti- 
mate can be made as to the accuracy of the approximation of (|6]) by comparing those 
correction terms with the HBT radii found by fitting correlation data with gaussians 
as in (|I^). If the former are much smaller than the latter, then (^) should be a good 
approximation. For example, the slopes and curvatures seen in heavy ion collision 
data generate 5Rf which typically have scales on the order of 

|ffi ^4(150MeV)2 ' (31) 

whereas R 2 , R 2 ± and R\ typically have scales on the order of R 2 ~ l/(75MeV) 2 |L4| 
so the approximations to these radii from eq. (|J) should be good to within roughly 
5%. As we will see later, however, the corrections could become important when 
determining the magnitude (and sign) of the cross term or the difference between R\ 
and R 2 S for systems with very short emission times. 

2.2 Boost Invariant Coordinates 

Now we would like to rederive the results of the preceding section using rapidities 
rather than longitudinal momenta, since the former boost invariant variables are 
usually more appropriate for relativistic collision experiments. Returning to eqn. (^), 
let us make an alternative on-shell definition of the 4-vector K\ 

K = (m t chY, K t , m t shF) (32) 

where K t = \{p u + p 2t ), m\ = m 2 + \K t \ 2 , Y = + y 2 ), and y* = |ln[(^ + 
PiL)/(Ei — PiLj]- Note that we use the subscript t throughout to denote transverse 
2-vectors as well as m t and other general transverse quantities; this should not be 
confused with the subscript _L which we use only to denote the "out" direction. 
Just as in the last subsection, we can expand the factor exp(iq-x) in eqn. (H) for 
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small momentum (and rapidity) differences. This time we find: 



C(y,q s ,q±,Y,K A 



1± 1-gW 



(33) 



q± [x —t ch(r] — Y)\ +jm t T sh(rj — Y) 

m t 



( <l_ ( x —t ch(r] — Y)j +ym t T sh(r] — Y) 



where y = yi— y 2 , r = \/t 2 — z 2 is longitudinal proper time, and// = ~ la[(t+z) /(t— z)\ 
is the space-time rapidity. (The reader should take care not to confuse the rapidity 
difference y with the cartesian coordinate y.) This time after exponentiating, we get 
a correlation function of the form: 



C(y, q s , q ± , Y, K ± )~l± exp -q 2 s R 2 s - qiR 2 ± - y 2 a 2 - 2q ± yR ±y 



(34) 



where again for the approximation of eqn. @ the correlation "radii" can be read off 
as the coefficients of the appropriate terms in eqn. (|34T): 



R 2 S 
R 2 , 



(y 2 ) 



a 



R 



J-y 



[x — (K±/m t )Tch(r] — Y)\ 
- (x — (K±/m t )Tch(r] — Y)) 2 
[m t Tsh(rj — Y)] 2 ^ — (m t Tsh(ri — Y)) 2 
{[m t x — K±Tch(r] — Y)] rsh(?7 — Y)) 
— (m t x — K±Tch(r] — Y)) (rsh(?7 — Y)) 



(35) 



Similarly to eqn. fl24|) , quadratic corrections which arise from expanding the 
denominator of @ for small q can be found to give 

Pi(pi)Pi(p 2 ) 



Jd*xS{x,K)\ 2 



exp ( -q 2 5R 2 - q 2 ± 5R 2 ± - y 2 5a 2 - 2q ± y5R ±J 



2f„,2 



(36) 
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where 



Y 



Note that these "side" and "out" corrections take the same form as those in eqn. 
(p5|), except that here rapidity rather than longitudinal momentum is held fixed 
while taking the derivative with respect to K±. Since experimental one particle 
spectra are usually presented as functions of rapidity and not longitudinal momentum, 
these corrections can be even more readily measured from the data than those of the 
previous subsection. 

Finally, we turn again to the corrections induced by using the exact off-shell 4- 
vector \{pi + P2) in the emission function of eqn. (0) rather than the approximate 
on-shell 4- vector K of eqn. (|32|). Making a Taylor expansion for small y and q t , we 
find 

\{Pi + P2) 

Furthermore, if we reparametrize the local 4-vectors of (|H]) in the following way 

u(x) = {u t ch£,ii p cos0, M p sin0, u t sh£) (39) 
where u t , u p and £ are independent of 0, then 



~ m 2 + \ql 1 - -4 + \q 2 s + \y 2 m 2 . (38) 



mt + hi I 1 



4) + k 2 + |y 2 

mi 4 4 



1/2 



u t ch(Y — £) — Kxjip cos ( 



(40) 
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Therefore, to quadratic order in y and q t 
Jd 4 xS[x,Up 1 +p 2 )}e i i" 



J d 4 xS(x, K)e i( i- X 



cxp 



n 1 -^? +9 - 2+y2m ' 2 



d 



dm 2 



In Pi (AT) 



(41) 



The most interesting thing to note about this off-shell correction is that it has no 
effect on the coefficient of the q±y cross term. 

Putting everything together, we find the following corrected model-independent 
expressions for the correlation radii of eqn. (|36|) : 

'lid Id' 



R Q 



(y 2 ) 



AK\ dK\ 2dm 2 



In Pi (it) 



P 



x — rch(?7 — Y) 

m t 



2 \ I \ 2 

) — { x rch(7? — Y) 

\ m t 



+ 



1 d 2 



K 2 S d 



4 dK'j_ 2 y m 2 dm 2 



lnPi(iT) 



a 



2 \ i 

[mtrsh(?7 — y)j ) — (mtTsh.(r] — Y))' 



+ I T 



1 C/ 2 1 , (j 



AdY 2 2 'dm 2 



lnPi(i?) 



R 



-Ly 



[m t x — K±Tch(r] — Y)] Tsh.(i] — Y) 



( mtX -K ± Tch(r ] -Y)){Tsh(ri-Y)) + ~^-^]nP 1 (K) (42) 

4 ar dK± 

Again, although R± y does not vanish in general, it does vanish for pairs with K± = 
(see appendix). 



15 



3 A Model with Nonrelativistic Expansion 

To get an idea of the usefulness of the model independent expressions just derived, 
we study a slight generalization of the thermal emission function presented in ||: 

S(x, K) = ^ exp H(t) I(p) J(z) . (43) 

Here T is a constant freeze-out temperature, and we define the space-time distribution 
of the source by a product of gaussians in the center of mass frame of an expanding 
fireball 

H(t)I(p)J(z) = , 1 exp - -4" - -^r I • (44) 



For the thermal smearing factor in eqn. (43), we take a nonrelativistic linear expansion 
4- velocity 



u{ x) 



I - (v R p/R G ) 2 - (v L z/L G ) 2 \ (1, v R x/R G , v R y/R G , v L z/L G ) , 
1 + \{vrp/Rg? + Uv l z/L g ) 2 , v r x/Rg, v R y/R G , v L z/L G ) , (45) 



where v R <C 1 and ^ C 1 are the transverse and longitudinal flow velocities of the 
fluid at p — R G and z = L G , respectively. 

Note that in the limit 5t — > 0, 5(x, K) becomes the Boltzmann approximation 
to the hydrodynamic emission function of eqn. (f|) with a constant freeze-out time to 
and a local chemical potential given by: 

/i(x) p 2 z 2 



T 2R 2 G 2L 2 G 



(46) 



In a sense, use of a nonzero 5t can be thought of as a smearing of the sharp 3- 
dimensional freeze-out hypersurface t = t over the fourth (temporal) dimension. 

Since the model is completely gaussian, analytic calculation of the one-particle 
distribution is straightforward, yielding 

W = (277)372 ex P ("T- + + ^ j ' (47) 
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500 1000 1500 
m t - m (MeV) 



2000 



Figure 1: The one particle spectrum of (47) is plotted as a function of m t — m for 
midrapidity (Y = 0) pions and kaons. The solid curve is for pions with no transverse 
flow (vr = 0), while the dashed curves are for pions (intercept normalized to 1) and 
kaons (intercept normalized to 2) with vr = 0.5c. The other source parameters used 
are R G = L G = 3 im, T = 150 MeV, and v L = 0. 
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Figure 2: The effect of longitudinal flow on the rapidity distribution of pions (outer 
two curves) and kaons (inner two curves) is shown for a source with no transverse flow 
vr = 0. The solid curves are for vl — while the dashed curves are for v = 0.5c. The 
other source parameters are defined as in fig. 1, and all curves have been normalized 
to 1 at Y = 0. 
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where 



1 

m 



R 2 G 



1 + 



E 



K 2 



R ' 



1 



T 2 



1 + 



E 



(48) 



In fig. 1 we plot P\ (K) as a function of m t — m for midrapidity (Y = 0) particles 
from a source with the parameters Rq = Lq = 3 fm, T = 150 MeV, and vl = 0. 
The decrease in the slope of the pion curves when the transverse flow is changed 
from vr — to Vr — 0.5c can be understood in terms of an effective blueshifted 
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T, 



eff 



1 + Vr 



From eq. 



it can also be seen that asymptotically as K\ 



oo 



Rl 



R 2 G T 



(49) 



(50) 



so that the K± dependence of the prefactor drops out and the spectrum takes the 
form of a pure exponential with an inverse slope of T e fj = 2T. Figure 1 also features 
a kaon distribution with flow which shows the same behavior. As can be seen from 
fig. 2, increasing the amount of longitudinal flow from vl = to Vl = 0.5c causes a 
widening of the rapidity distribution both for pions and kaons. 
Using (|3"0"|) the correlation radii are readily found to be: 
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where the corrections due to using eqn. (0) rather than eqn. (Q) have been grouped in 
the square brackets, and we have only kept terms up to second order in the velocities 
in order to be consistent with the nonrelativistic approximation ([45]). 

First we would like to note that by rotating the coordinate system for each pair 
to z = K (so that K± = by definition) and neglecting the correction terms, our 
expressions for the correlation radii reduce to those of || for the nonrelativistic (E^ = 
m) and spherically symmetric case of Rc = L G and vr/Rq = vl/L g = l/t . As they 
point out and can be seen from eq. (f48[), transverse flow causes the "side" radius to 



measure something smaller than the real geometrical radius Rq- In fact since this 
is a completely gaussian model, it should be no surprise that R s simply measures 
the length of homogeneity of eq. fll~9|) in the transverse direction. Similarly, Li just 
measures the longitudinal length of homogeneity which will be smaller than Lq if 



vl > 0. One of the most interesting features of ([51]) is that in the absence of the 
corrections, not only is the difference between R\ and R 2 S directly proportional to 
the square of the emission time (which is simply the length of homogeneity in the 
temporal direction), R\ L is as well. 

Notice that for systems with T ~ 150 MeV and R s , R±, Rl ~ 3 fm, the correction 
terms do not alter the naive expressions for those radii by more than 3%. For very 
small radii and very short emission times St, however, the correction terms may 
actually have a noticeable cancellation effect both on the magnitude of the cross 
term and on the difference between R\ and R 2 S . This should be kept in mind when 



extracting limits on St from the data |14[ . For example, for pions with K± ~ m, this 
kind of cancellation will occur for emission times St < \ fm. In particular, for St = 0, 
R\ would actually be smaller than R 2 S in this model. However, since present heavy 



ion correlation radii are measured to be around 3 fm |14| and the experiments are 
not yet able to resolve 3% effects, keeping the correction terms may not be necessary 
when comparing a specific model to heavy ion correlation data. 

One might at first think that the cross term for this model would vanish if the 
radii were calculated in the LCMS frame, since (3l = in that frame. This is not 
the case, however, because the emission function S(x, K) is not longitudinally boost 
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invariant, even in the case of non-relativistic Galilei-transformations. After making 
the appropriate transformations into the LCMS frame 



t' = i L {t- f3 L z) z' = lL (z - [3 L t) 



1L 



-1/2 



(52) 



t' z' cross terms are introduced into the gaussians. These in turn give rise not only to 
a nonzero R\ L cross term but also modifications to the other radii. Neglecting the 
correction terms, 
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(53) 



where /3l and 7^ in the above expressions are evaluated in the fixed center of mass 
frame. Note that in this frame there is now also a geometrical contribution ~ L 2 to 
R' 2 L and R'l — Rf; since it is multiplied by a factor 1/c 2 relative to the (St) 2 terms, it 
vanishes in the non-relativistic limit c — > 00. However, the (St) 2 contribution to R'l L 
in particular survives in this limit. 



4 A Model with Relativistic Longitudinal Expan- 
sion 

Now we move to a model similar to those in [|K| which should provide a more realistic 
description of particle emission from a relativistic collision. In the center of mass frame 
of an expanding fireball, we define the following emission function 
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(54) 



where again T is a constant freeze-out temperature, r = \/t 2 — z 2 is the longitudinal 
proper time, and Y — \ \v^E K + K L )/(E K — K L )] is the rapidity of a particle with 
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momentum K. This time in the limit 5r — > 0, (54) becomes the Boltzmann approxi- 
mation to (f|) with a constant freezeout proper time r and a local chemical potential 
given by 

Ms) _ p 2 ^ 

T " 2i% 2(^) 2 1 J 

The second exponential in the emission function (|54j) can be interpreted as the 
space-time distribution of point-like sources, each of which emits a thermal spec- 
trum, boosted by the flow 4- velocity u(x), as given by the first exponential and the 
ch(i] — Y) prefactor. For simplicity, the source distribution in space-time is taken to 
be gaussian. 

For this model, we consider a flow which is still non-relativistic transversally but which 
now exhibits Bjorken expansion (fluid rapidity = space-time rapidity) longitudinally, 

u(x) ~ ((l + \{vp/R G ) 2 ) chry, (vx/R G ), (vy/R G ), (l + \{vp/R G ) 2 ) sh V ) , (56) 

where v <C 1 is the transverse flow velocity of the fluid at p = R G . This flow profile 
corresponds to a longitudinal velocity VL(z,t) — z / 1 . With this definition, K-u takes 
the following longitudinally boost-invariant form 

K-u = m t [l + ^vp/Rg) 2 }^ -Y)- K ± (vx/R g ) . (57) 

If we restrict ourselves to particle pairs with m t >T and \Y\ <C 1 + (Sr]) 2 mt/T, 
then we can perform a modified saddle point approximation by expanding ch(r] — Y) 
in (|57j) in powers of rj' = rj — Y, keeping in the exponent only terms up to second 
order and expanding everything else to the desired order. For our calculations, we 
approximate S(x, K) by 
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Note that we keep only the r]' 2 term when expanding the chr/ prefactor, but we also 
keep a term {m t /T)r]' A from the expansion of the exponent. The latter term is to be 
taken as be roughly of the same order as r/ 2 for reasons which will become clear later. 

For the particle emission time, it must be true physically that 5t/tq < 1. Rather 
than demanding the much stricter condition 5t/t$ <C 1, we simply assume that 
this ratio is small enough (e.g. ^t/tq<|) so that we can replace integrals over 
only positive values of r with ones ranging from — oo to +00. Finally, in all of our 
calculations we throw away all terms of 0(f 4 ), in keeping with our nonrelativistic 
approximation in the transverse direction. 

Given these approximations, calculation of the one particle distribution can now 
be done analytically, yielding 
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and our expansion parameter is defined by 
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Note that for pairs in which m t /T ^> l/(Sr)) 2 as were studied in ||20|| , {8rj)1 becomes 
simply T/m t . This is the reason that we consider (m t /T)(5r])l to be of the same 
order as (5rj) 2 . 

In fig. 3 we plot numerical calculations of P\{K) as a function of m t — m for 
midrapidity (Y = 0) particles from a source (|54f) with the parameters To = 4 fm/c, 
Rg = 3 fm, 5i] = 1.5, and T = 150 MeV. We have checked that our analytic expres- 
sions provide excellent (< 5% error) approximations to the exact numerical results. 
Again the decrease of the pion slope as the transverse flow parameter v is increased 
from v — to v — 0.5c can be well understood in terms of the effective blueshifted 
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Figure 3: The one particle spectrum obtained by numerically integrating (54) is 
plotted as a function of m t — m for midrapidity (Y = 0) pions and kaons. The solid 
curve is for pions with no transverse flow (v = 0), while the dashed curves are for 
pions (normalized to 1) and kaons (normalized to 2) with v = 0.5c. The other source 
parameters used are r = 4 fm/c, Rq = 3 fm, 5i] = 1.5, and T = 150 MeV. 
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Figure 4: The rapidity distribution of pions (solid curves) and kaons (dashed curves) 
is shown for a source with no transverse flow v = and other source parameters equal 
to those of fig. 3. 
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temperature of eq. (|49|). In contrast to the nonrelativistic case, however, Pi does 
not quite become a pure exponential as K± — » oo, because although the prefactor Rl 
again cancels the prefactor m t , there is an additional K± dependence coming from 
the prefactor 



(<fy), -> ^T/m t . (62) 

The effect of this prefactor can also be seen at the origin, where its partial cancellation 
of the m t prefactor results in less curvature than is seen in the nonrelativistic case of 
fig. 1. For completeness, in fig. 4 we also show the rapidity spectrum both for pions 
and kaons. Although there is a slight decrease in the width for kaons, the effect is 
much smaller in this case than in fig. 2 due to the relativistic longitudinal flow which 
causes the difference in mass to become less important. 

The correlation radii can now be calculated by using eqs. ( |5UD or P2]). By means 
of a saddle point approximation, the correlation radii are expanded to the second 
order in the small parameters (5rf)l and (St/tq) 2 . Therefore, in performing these 
calculations below, we only keep terms of order (Srj)* (or (m t /T)(5r))®) except for 
smaller terms involving (5t/t ) 2 or v 2 , in which case we keep only terms of order 
(5r))l and (Srj)®, respectively. Please note that, since the general expressions (|30|) or 
( fl"2~D for the correlation radii are all at most quadratic in the proper time, no term 
involving (<5t/t ) 4 will in fact occur in the results. 

4.1 HBT Radii in Cartesian Coordinates 

In cartesian coordinates, the correlation radii take the following form (ordered by 
powers of the small expansion parameters (5r/r ) 2 and (5rj)l): 
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(63) 



Here v — 1 + (R*/Ra) 2 — ^( m t/T)(5f])l, and we have neglected the corrections which 
come from using (|2]) instead of (||). Although this model is not completely gaussian, 
within the scope of our approximation R 2 still roughly measures the transverse region 
of homogeneity of the fluid, as can be seen by comparing fl60|) with (19). Although we 
will show that in practice all terms given in (|63| ) are important, we will for didactical 
purposes first consider only the leading order in the small expansion parameters. 
Then the expressions (|63|) simplify and can be reformulated as follows: 
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K x shY 

In agreement with [JTO]] we find that, in each principal direction of the expanding fire- 
ball, two length scales should be distinguished. In addition to the geometric length 
scales Rq and Lq = r 5r] in the transverse and longitudinal directions respectively, 
we have two "lengths of homogeneity" generated by the flow gradients. The transver- 
sal and longitudinal homogeneity lengths are given by the following expressions: 
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Please note that the occurrence of Tq in both longitudinal lengths, Lq and Ljj, has 
two different origins: whereas the geometrical longitudinal extension of the fireball 
at freeze-out is clearly always proportional to the mean freeze-out proper time To, 
its occurrence in the longitudinal homogeneity length is due to the specific choice 
of the velocity profile, since for a longitudinally boost invariant velocity profile the 
velocity gradient is just given by the inverse proper time. In fact, the true origin of 
the homogeneity lengths ( p5f) is seen by writing them in the form 

T 1 
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H m t (dvt/dp) 2 ' 

(66) 
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where Vt(p) = vp/Rc, and in the second line = (chrj, 0, 0, sir?/) denotes the lon- 
gitudinal part of the flow velocity profile (|56|) which satisfies d^u^ — 1/r. Eq. ( |66| 
makes the nature of the homogeneity lengths explicit in showing how they are gener- 
ated by the flow gradients at freeze-out. 

With these notations, the correlation radii can be written as follows : 
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As already pointed out in |T(J, the correlation radii are seen to be dominated by the 
shorter of the geometric and homogeneity lengths. This means in particular that if 
v 7^ 0, then R 2 S will be smaller than the geometrical radius Rq- As the transverse 
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mass increases, this reduction of R s and R± relatively to the pure geometric radius 
becomes more pronounced. In the longitudinal direction, as a general consequence 
of the particle pair motion with velocity Y, the system appears Lorentz-contracted. 
Hence g^-correlation functions at finite values of Y measure longitudinal correlation 
radii, which are reduced by the corresponding Lorentz-contraction factor ch _1 V, as 
shown by (p7|). Similar purely kinematic factors affect the out-longitudinal radius 
R±l- 

Returning now to the higher order corrections shown in (p3| ) we observe that, in 
contrast to the nonrelativistic model, the difference between the squares of the "out" 
and "side" radii depends on the rapidity Y (or fa) of the pair and is not quite directly 
proportional to the duration of particle emission (St) 2 even for pairs with (3l = ||21|| . 
It is also worth noting that although R\ L vanishes when K± = or Y = 0, for high 
K± and \Y\ pions R\ L is of the same order of magnitude as R\ y so it has a significant 
effect on the form of the correlation function. 

This can be seen most easily in a numerical example. For simplicity, we consider a 
pion source with no transverse flow (v = 0) which freezes out instantaneously (St = 0) 
with the following other source parameters: Rq = 3 fm, tq = 4 fm/c, Sr] = 1.5, and 
T = 150 MeV. Given these parameters and any set of momenta q, K, it is possible to 
determine the correlation function both by using the approximate radii of ( p3l) and 
by performing an exact numerical calculation of the correlation using (0) and (J54[) . In 
all of our plots of the correlation function, solid curves are used to denote numerical 
calculations, while dashed curves are used to denote our analytic approximation. 

The symmetric curves in fig. 5 show the correlation as a function of qi for Y = —2, 
q± = 30 MeV, q s = and K± = 0. In this case as we mentioned earlier, since K± = 
the cross term vanishes and the correlation function peaks at qi = 0. As K± is 
allowed to increase, however, R\ L causes the peak to shift toward negative values of 
qi, as can be seen in the asymmetric curves which have been calculated for K± = 200 
MeV and all of the other momenta the same. Similarly, if Y is allowed to increase 
to 0, the maximum shifts back to qi = 0, and for Y > 0, the maximum is located 
at a qi > 0. It should also be pointed out that in each of the above cases, the 
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Figure 5: The correlation function generated by the pion source of (54) is plotted as a 
function of qL for two values of K±: K±_ = (symmetric curves) and K± = 200 MeV 
(asymmetric curves). The solid lines are exact numerical results, while the dashed 
lines are our analytic approximation (15,63). The source parameters are v — 8r — 0, 
Rg = 3 fm, To = 4 fm, 8r] = 1.5, T = 150 MeV. The pair momenta we fixed at 
Y = -2, q ± = 30 MeV and q s = 0. 
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direction of the shift of the peak is reversed if a negative q± is used instead of a 
q± > 0. To give a quantitative idea of how good the analytic approximation is for 
the K± = 200 MeV case, we found that the best gaussian fit to the numerical curve 
could be reproduced by multiplying the correlation radii by the following factors: 
R 2 ± -> 0.92R 2 ± , Rl -> 0.95R 2 L , and R\ L -> 0.75^. 

Figure 6 shows the correlation as a function of qj_ for two different values of q^. 
Both sets of curves are calculated for Y = —2, q s = and K± = 200 MeV, but the 
upper ones have qi, = while the lower (asymmetric) ones are for qi = 100 MeV. It 
can be seen that increasing q^ from has the effect of shifting the peak down and to 
the left (to q± < 0). This figure shows clearly that interesting physics could be missed 
if correlation models are only plotted as a function of a single momentum difference 
with all other set equal to zero. Again to get a quantitative idea of the validity of 
the analytic approximation, we found that the best gaussian fit to the numerical curve 
for q L = 100 MeV could be obtained using the factors R\ -> 0.85^1, R\ -> 1.08-Rf,, 
and R\ L -> 0.75R 2 ±L . 

As can be seen from figs. 5 and 6, the simple analytic expressions of fl63|) repro- 
duce the exact correlation functions remarkably well considering the crudity of the 
approximation. By extensively exploring the parameter space of the model, we have 
found that the quantitative error estimates we have obtained in figs. 5 and 6 are 
somewhat typical of the maximum discrepancies for reasonable parameters. Namely, 
the analytic approximations of (|63]) for Rj_, R\ and R\ L are able to reproduce the 
best gaussian fits to the numerical expressions to within <20%, <10%, <33%, re- 
spectively (e.g. for Rj_ L , (1 — .75)/. 75 ~ 33%). Although not shown, the analytic 
expressions for R 2 S are much better, their discrepancy from numerical fits is typically 
<5%. We would also like to note that we have performed numerical calculations using 
eq. (H) and find them to agree to within 3% with numerical calculations using (0), so 
we are well justified in neglecting those corrections in eqs. fl63|). 

The analytic expressions of fl63| ) are even better approximations for heavier parti- 
cles like kaons, since for them m t /T > 3 so (Sr])l forms a smaller expansion parameter. 
This behavior can be seen in fig. 7 where we plot the kaon correlation as a function 
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Figure 6: The same source parameters as in fig. 5, but with pair momenta Y = —2, 
K± = 200 MeV, q s = 0, are used to plot the correlation as a function of q± for q L = 
(symmetric curves) and qi = 100 MeV (asymmetric curves) 
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Figure 7: The narrower curves show the kaon correlation as a function of q± from a 
source with the same parameters as in fig. 5 but with momenta defined by Y = +2, 
K± = 200 MeV, q s = and qi = 150 MeV. The wider curves have been obtained by 
using v — 0.5 instead of v — 0. 
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of q± for Y = +2, K± = 200 MeV, q s = and = 150 MeV. The narrower curves 
are obtained by using the same source parameters as in figs. 5 and 6, while the wider 
curves feature a transverse flow parametrized by v — 0.5c. The best gaussian fit 
to the wider numerical curve can be obtained in this case by multiplying the radii 
of the wider analytical curves by the factors R\ — > 0.94i?^_, R\ — > 1. 02-Rf,, and 
Rl L - 0.9R 2 ±L . 

Perhaps the best way to study the correlation function is to make a 2-dimensional 
surface plot of C — 1 as function of both q± and qi- Figure 8 shows such a plot of 
the numerical calculation of C — 1 for Y = —2, K± = 200 MeV and q s = 0. The 
effect of the cross term can be seen in the form of a ridge running from the peak 
at q± = qi = down to the front left where qi > and q± < 0. Since cylindrical 
symmetry precludes the existence of "side-out" or "side-longitudinal" cross terms, the 
only effect of averaging over q s from to some maximum value such as 30 or 50 MeV 
HU would be to reduce the intercept of the correlation function to some value less 



than 1. This averaging, however, should have very little impact on the qualitative 
ridge structure of the "out-long" correlation function. Consequently, this kind of ridge 
should be clearly identifiable experimentally and in fact may have already been seen 
in preliminary E802 correlation data G^[. 

Before analyzing this model in rapidity coordinates, we would like to note that 
the LCMS radii of this model can be obtained simply by setting (3l = and Ek = vn t 
in (|63"D . Note that the factor of Y in R\ L should not be set equal to zero, since it 
arises from the space-time rapidity distribution of the point-like sources in ([54]) which 
obviously breaks the boost invariance of the emission function in the longitudinal 
direction ||23|| . Transforming to the LCMS frame introduces a Y dependence which 
eventually translates into a nonvanishing cross term. We would like to emphasize 
that, to the first order in the small expansion parameters, our results reduce to the 



expressions for the LCMS correlation radii derived in [10]. However, in the light of 
a comparison of the results obtained within the framework of our analytical approx- 
imation with an exact numerical computation of the correlation function, it turns 
out that the second order contributions to the correlation radii must be included. In 
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Figure 8: The numerically calculated correlation function generated by the pion 
source of fig. 5 is plotted as a function of qi and q± for Y = —2, K± = 200 MeV, 
and q s = 0. 

particular, the out-longitudinal cross terms, whose effect can be clearly seen in Fig.8, 
is completely missed at leading order. Nevertheless, for this model, the LCMS frame 
has the advantage that the expressions for the correlation radii are much simpler than 
for those in the fixed frame. On the other hand, this same simplicity can be achieved 
without the complication of reference frame shifting by expressing everything in terms 
of boost-invariant coordinates, as we will now show. 

4.2 HBT Radii in Boost-Invariant Coordinates 



Using the model independent expressions of ([42") along with the emission function of 
), we obtain the following correlation radii 
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where in contrast to section 4.1, Y is now defined Y — |(yi + y 2 )- Note also that in 
contrast to the corresponding radii of (|63D R± and a in the above approximation are 
both independent of rapidity. In addition, the cross term R± y will be small compared 
to these radii, especially for higher mass particles like kaons which have (Sr])l •€ 1 or 
for future ultrarelativistic collisions in which (Srj) 3> 1. 

The astute reader will note that aside from a difference in the definition of Y, the 
fixed frame correlation radii of the last subsection can be easily derived from those 
of flBq) in the following way: First insert the radii of (|68|) into the expression (|34|) 
for the correlation function, then make the replacement y — > qi/Ex — (3LK±q±/mf, 
rewrite the resulting expression in the form of eq. fl!5D , and finally read off the radii 
of eq. fl63|) . The reason for this can easily be seen by noting that 

q-x ~ q± — — r ch(r? — Y) — y m t r sh(r? — Y) — q±x — q s y 

m t 

~ q±—r ch(?7 -Y)- (^qL ~ —PLq±) t sh(r/ - Y) - q ± x - q s y (69) 
m t \E K m t J 

where in the top line Y = |(yi + y 2 ), while in the bottom line Y = \\o[{Ek + 
Kl)/(Ek — Kl)}- Note that in particular the LCMS radii can be found simply by 
making the replacement y — > qL/m t . Based on this equivalence, one can see that 
for systems undergoing Bjorken longitudinal expansion, LCMS correlation functions 
are nothing more than approximations of fixed frame correlation functions in rapidity 
coordinates. Since the latter formulation is manifestly boost invariant and avoids the 
complications arising from the introduction of the different LCMS-reference frames, 
it is much more desirable to use those coordinates. For the remainder of this section, 
we use the definition Y — |(yi + y 2 )- 

Using the same source parameters as in the last section, figure 9 shows the pion 
correlation as a function of y for Y = —2, K± = 200 MeV and q s = 0. The sym- 
metric curves are for q± = while the asymmetric curves are for q± = 30 MeV. The 
best gaussian fit to the asymmetric numerical curve can be obtained in this case by 
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Figure 9: The same source parameters as in fig. 5 are used to plot the correlation as a 
function of y = y x — y 2 for q± = (symmetric curves) and q±_ = 30 MeV (asymmetric 
curves). For both curves the pair momenta have been fixed to be Y = —2, K± = 200 
MeV, and q s = 0. 
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multiplying the radii of the wider analytical curves by the factors R\ — > 1.17i?5_, 
a 2 — > 0.96a 2 , and i?^ y — > 1.6i?j_ y . Again, these corrections are somewhat typical of 
the maximum discrepancies, and the analytic approximations of for Rj_, R 2 , a 2 
and R± y are thus able to reproduce the best gaussian fits to the numerical expressions 
to within <20%, <5%, <10%, <40%, respectively (e.g. for R ±y , 0.6/1.6 ~ 38%). Of 
course, much better agreement could be obtained if a more sophisticated analytical 
approximation is used in place of eq. (|58|). 

5 Conclusions 

By taking second derivatives of the two particle correlation function around q = 0, we 
have derived model-independent expressions for correlation radii both in cartesian and 
boost-invariant momentum coordinates. In both cases, an "out-longitudinal" cross 
term arises naturally. In the context of two "gaussian" models, this term is found 
to have a significant effect on the form of the correlation function. We therefore feel 
that future correlation data should be fit to one of the following two functions 

C(pi, p 2 ) = 1 ± A exp (-q 2 R 2 s - q\R\ - q 2 L R 2 L - 2q ± q L R 2 ±L ) (70) 

or even better 

C(p u p 2 ) = 1 ± A exp (-q 2 s R 2 s - q 2 ± R 2 ± - yV - 2q ± yR ±y ) , (71) 

where R\ L (or R± Y ) can be either positive or negative. 

Currently, data is usually fit to ([F(]) with R\ L a priori set equal to zero [14-17]. 
For that reason, all qi with the same \qi\ are usually binned together, but in the 
process, the relative sign between q± and qi gets lost. This procedure effectively 
averages out the cross term at the expense of introducing large systematic errors into 
the measured "out" and "longitudinal" radii. One way to avoid this averaging in 
practice is to define the ordering of particles 1 and 2 by always demanding that q^ (or 
y) be positive. This then determines the sign of q± (and q s ) so that positive values 
can be binned separately from negative values, allowing one to generate plots like the 
one we have shown in fig. 8. Not only will measurement of cross terms provide new 
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information about the emitting source, it should greatly increase the accuracy of R± 
and Rl (or a) measurements. 

The mo del- independent expressions for the radius parameters of the HBT corre- 
lation function show very clearly that these parameters do not generally measure the 
geometric size of the source, but rather its lengths of homogeneity in the four space- 
time directions. For expanding sources like those created in heavy ion collisions, the 
gradients of the thermodynamic parameters and of the flow velocity field contribute 
to the inhomogeneity of the source. In fact, regions of homogeneity may extend over 
only a small fraction of the source, in which case the two-particle correlation function 
is sensitive only to these subdomains. Moreover, particle pairs with different average 
momenta will generally see regions of homogeneity with different size, giving rise to 
a characteristic independence of the correlation radii. 

In this paper we have studied these features quantitatively for sets of cylindrically 
symmetric models with gaussian density profiles in which the sources undergo longi- 
tudinal and transverse collective expansion but freeze out at a constant temperature. 
The effect of the flow gradients on the lengths of homogeneity and on the spatial HBT 
size parameters has been seen explicitly. They lead to a reduction of the correlation 
radii relative to the geometric radius parameters, and this effect increases with the 
average momentum of the pair relative to the center-of-mass of the source. The tem- 
poral length of homogeneity of the source, given by the duration 5t of the emission 
process, affects both the difference R\ — R 2 S (as has been noted previously @, [J) 
and the new "out-longitudinal" cross term. The effects of possible gradients of the 
freeze-out temperature have not yet been studied in this context, but are expected 
to have similar qualitative consequences. In fact, a difficulty in separating effects 
of flow gradients from those of thermal gradients was noted before in the context 



of a spherically symmetric model [19. El. It was found that both mechanisms can 



lead to a concave curvature of the single particle m r spectra [IS], as well as a similar 
K±_- dependence of the "side" and "out" radii in the HBT correlation function p3] . 

In [Rl |10| the difference between the geometrical and HBT radii has been expressed 
in terms of a so-called "thermal radius" . Our analysis shows that it is really not the 
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existence of a temperature, but of a flow velocity gradient which causes the appearance 
of a length of homogeneity in the HBT radii. The temperature only plays a role as a 
smearing factor, and the ratio T/mt sets the scale at which the inhomogeneity of the 
flow field becomes effective. Different flow velocities in the transverse and longitudinal 
directions generally lead to different transverse and longitudinal homogeneity lengths, 
Rh and L H . In || ffO this was not obvious because the flow gradient was fixed to be 
1/tq in all directions by the choice of the flow velocity profile. 

All of our calculations in this paper were done in a fixed reference frame, thus 
avoiding the complications with the LCMS frame discussed in Section 2. However, 
we found that by parametrizing the correlation function in terms of rapidities rather 
than longitudinal momenta, one finds a longitudinal correlation radius and an out- 
longitudinal cross term which for sources with boost-invariant longitudinal expansion 
can be well approximated by the LCMS results. Since this parametrization avoids 
the LCMS problems of shifting frames, we suggest that the concept of the LCMS be 
abandoned in favor of using rapidity coordinates. We also showed that the existence of 
an out-longitudinal cross term is not affected by this choice of coordinates or frames, 
although its actual size is. 

The analytic expressions for the HBT size parameters developed in this paper have 
been tested numerically and were found to be sufficiently accurate for being useful in 
obtaining good qualitative insights on the effects which various features of the source 
have on the shape of the correlation function. We also studied explicitly the usually 
neglected corrections due to the off-shell nature of the average 4-momentum entering 
in the correlation function and found them to be very small (< 3%). To the extent 
that our two models for the source emission function are reasonable approximations to 
reality, these relations can be used to study the effects of longitudinal and transverse 
flow and of the time and duration of the freeze-out process on the HBT data. We 
have checked that the models produce single particle spectra with reasonable shapes 
which very likely can be used for good fits to the data (in particular once resonance 
decays are included). A more detailed analysis of the HBT data in the framework of 
these models thus appears as an attractive project. 
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6 Appendix 



In this appendix, we prove that in the limit K± — > 0, R 2 ± — > R 2 and the cross term of 
either eq. (0) or eq. vanishes (depending on the coordinate system used). The 
crucial ingredient of the proof is that the emission function is a Lorentz scalar whose 
K dependence only enters in the form of scalar products with cylindrically symmetric 
local 4-vectors. For simplicity in the following, we will assume that there is only one 
such local 4-vector, but generalization to the emission function of fliTf ) can be done 
trivially. 

We assume an emission function of the form 

S(x,K) = S(t,p,z } m 2 } iP) (72) 

where in a cartesian coordinate system 

ip = K-u{x) = Ek uo — K± Up cos <p — Kl u z (73) 

and uo, u p , and u z are independent of (p. Using rapidity coordinates as in (j3^) and 
59f) on the other hand, we can see that 



if) = m t u t ch(Y — £) — K± u p cos(fi (74) 

where u t , £ and u p are independent of 4>. 

In either long as the ip dependence of S is smooth, it follows that 



lim^ J d A x coscfi f(t, p, z) S = 



f OS 
lim / d A x cos 6 fit, p. z) — — = 
k^oJ rM,n 'dip 



r a S 

lim / d 4 x cos 6 f (t, p. z) — — - = (75) 
k^oJ dip 2 

for any ^-independent function /. From the first of the above equations we can see 
that 

lim(/(t,p,z)cos0) = O (76) 
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in particular (x) = (pcoscp) = 0, and consequently the non- derivative terms of R\ L 
and R± y vanish in the limit K± — > 0. Furthermore, since 

lim (p 2 cos 2 <P) = lim (p 2 sin 2 0) = Up 2 ) , (77) 

i.e. (x 2 ) = (y 2 ) = ^(p 2 ), it can be seen that the non-derivative terms of R\ equal 
those of R 2 in that limit. 

As for the momentum derivative terms, from eq. ( |7o] ) we have 

d r (99 

Jro^w--/^^ 008 ^" (78) 

in either set of coordinates. Similarly, 

k ± ^o dK L dK± k ± ^o dY dK± 

Since 

JLjU^*, = ^^^y?) _ ( dP ^\ (" p ^ 



dK L dK ± P\{K) dK L dK A 



Pi{K) 



dKr dK\ 



(80) 

we have proved that R\l of fl30|) vanishes for K± — > 0. The proof for R± y follows 
simply by replacing d/dKi with d/dY in the above equation. 

The momentum derivative term for the "out" radius does not vanish in this limit, 
rather in the cartesian system it takes the form 

d 2 „. r> , r . ( u dS 2 2A d 2 S\ 



lim ^{11)= lim ^M^ + ^cos 2 , 
K ± ^odK 2 k ± ^oJ \Ek dip p 



E K dip p T dip 2 J 

The derivative term for the "side" radius is a bit trickier since it involves a ratio of 
two quantities which vanish in the K± — > limit 

1 ^ i t>\ l- _ f j4 ( u o u p j\ 9S 



lim — — — PAK) = lim / d 4 x - cos0 — (82) 
K ± ^oK ± dK ± n ; k x -*oJ \E K K± V dip y ' 

Determination of the appropriate limit of the second term above is found by the rule 
of l'Hospital by dividing the derivative (with respect to K±) of the numerator by the 
derivative of the denominator. When this is done, the results for the "side" and "out" 
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directions become identical. A similar argument can be used to show the same thing 
in the rapidity coordinate system. 
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